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FLUCTUATIONS OF THE EULER-POINCARE CHARACTERISTIC 
FOR RANDOM SPHERICAL HARMONICS 

V. CAMMAROTA, D. MARINUCCI AND I. WIGMAN 


Abstract. In this short note, we build upon recent results from [7j to present a precise expression for the 
asymptotic variance of the Euler-Poincare characteristic for the excursion sets of Gaussian eigenfunctions 
on S 2 . 


1. Introduction and main result 

The geometry of excursion sets for Gaussian random fields has been a subject of intense research over 
the last fifteen years; much work has focussed on the investigation of the Euler-Poincare characteristic, 
henceforth EPC [3]. We recall here that the EPC %(A) is the unique integer-valued functional, defined on 
the ring C of closed convex sets in R w , which equals x(A) = 0 if A = 0, x(A) = 1 if A is homotopic to the 
unit ball, and satisfies the additivity property 

X{A U B) = x(A) + x(B) - x(A n B), for all A,B eC. 

Clearly, the EPC is a topological invariant (i.e. it is invariant under homeomorphisms); its investigation for 
the excursion sets of random fields was initiated in the late seventies by Robert Adler and his co-authors. 
This stream of research has eventually resulted with the discovery of the beautiful Gaussian Kinematic 
Formula (GKF) pHE]. 

More precisely, let / be a real valued random field defined on the parameter space JA\ its excursion sets 
are defined as 

A u (f-, M) = {x € M : f(x) > u} , u G R. 

Let Cj ’s for j = 0,..., dim(Ad), denotes the Lipschitz-Killing curvatures for the manifold Ai with Riemann- 
ian metric g* induced by the covariance of /, i.e., for U X ,V X £ T X M, the tangent space to A! at a; we 
have 

gl(U x , v x ) := E[(C4/) • 04/)], 

(see [2] for further details); in particular Co is the EPC. The functions pj ’s are the so-called Gaussian 
Minkowski functionals and they are defined by 

(1.1) pj(u) = (2tt )-^/2H 3 _ l{ u)e- u2 /\ 
where H q (-) are the Hermite polynomial of order q : 

dp 

H-i(u) = 1 - $(u), Hj(u) = (-1)- J (0(u)) -1 ^— (j)(u), j = 0,1,..., 

</>(•), $(•) denote the standard Gaussian density and distribution functions, respectively; for example: 

Ho(u) = 1, Hi(u) = u, H 2 (u) = u 2 — 1, Ho(u) = u 3 — 3 u ,... 

The GKF states that the expected EPC of the excursion sets of a smooth, centred, unit variance, Gaussian 
random fields / : A4 —> R is 

dim(M) 

(1.2) E [ X {A u {f-M))\= E ^(M) Pj {u). 

3=0 

While the GKF yields a precise expression for the expected value of the EPC of excursion sets of smooth 
Gaussian processes, the analysis of higher moments, and, in particular, of the variance, is still open. The 
latter question is of both theoretical and applied interest; for instance, in the recent paper |T) five different 
methods are suggested to estimate numerically the covariance matrix of the EPC characteristic for the joint 
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excursion sets at various thresholds. These results were subsequently exploited to approximate excursion 
probabilities, the so-called Euler-Poincare heuristic (3] Section 5.1. 

In this paper, we establish analytic formulae for the covariance of the EPC characteristic of excursion sets 
at different thresholds, focussing on an important class of fields: Gaussian spherical harmonics. We establish 
a rather simple expression which seems to be closely related to a second-order Gaussian Kinematic formula, 
in a sense to be made clear below. More precisely, consider the Laplace equation 

A s 2 ft ~ = 0) ft '■ S 2 ~t 

where A 52 is the Laplace-Beltrami operator on S 2 and \i = —£(£ + 1), £ = 0,1, 2,.... For a given eigenvalue 
X(, the corresponding eigenspace is the ( 2 £+ l)-dimensional space of spherical harmonics of degree £; we can 
choose an arbitrary L 2 -orthonormal basis {Yi m (.)} m __ e and consider random eigenfunctions of the form 

1 1 

h{x) = / 2£ + 1 a trn Y tm(x), 


where the coefficients {ai m } are independent, standard Gaussian variables. The law of fi is invariant w.r.t. 
the choice of a L 2 -orthonormal basis {F; ra }. The random fields {fi(x), x £ <S 2 } are centred, Gaussian and 
isotropic, meaning that the probability laws of fi(-) and fi(g-) are the same for any rotation g £ SO( 3). 
From the addition theorem for spherical harmonics m Theorem 9.6.3) the covariance function is given by 


E[fi{x:)fe(y)] = P e (cosd{x,y)), 

where Pi are the Legendre polynomials and d(x, y) is the spherical geodesic distance between x and y. An 
application of the GKF (1.2 1 gives in these circumstances: 

£{£ - 


(1.3) 


E[; x{A u (fi;S 2 ))} = ^Eexp{-u 2 /2}u- 
a/tt 


1 ) 


2[1 -$(«)], 


for a proof of formula (1.3) see, for example, uni Lemma 3.5 or m Corollary 5. Note that, as u —> — 00 , 
the right hand side of (1.3) yields the Euler-Poincare characteristic of the two-dimensional sphere i.e. 

lim E[ x (A u (U,S 2 ))} = 2. 


The analysis of spherical Gaussian eigenfunctions is motivated by applications arising mainly from Math¬ 
ematical Physics and Cosmology. In particular, Gaussian eigenfunction have been conjectured [0 to ap¬ 
proximate deterministic eigenfunctions on generic billiards (surfaces with smooth boundaries). On the other 
hand, spherical Gaussian eigenfunctions are the Fourier components of isotropic spherical random fields, 
and, because of this, have been deeply exploited in the analysis of cosmological data, see for instance [ 12 ]. 
Let I C E be any interval in the real line and 

Ai(ff,S 2 ) = fi\l) = {x G S 2 : ft(x) £ 1} . 

Our principal result is the following: 


Theorem 1. As £ —> 00 , for every intervals I\,l 2 Q R, 

(1-4) Cov[x(4(//;5 2 )), x(A l2 (f£,S 2 ))} = 1 ^X 2 + 0(£ 5 / 2 ), 

where Ifor i = 1,2, are given by 

Pi = J p(ti)dti , p(t) = (- 1 4 + 4t 2 - l)e“T. 

The constant involved in the O(-) notation is universal. 


In particular the variance for any given interval follows as an easy corollary: 


Corollary 1. For every interval I CB, 

(1.5) Var [ X (Mf ( -,S 2 ))\ = ^ 1 2 + 0{P' 2 ) 

where 

I = J (—t 4 + 4 1 2 — l)e - ^" dt. 
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(a) positive range 


(b) negative range 


Figure 1. A_(-if + 4tf - + 4t| - l)e _ ^e“ 


As explained in Section[3j the proof of Theorem [l] follows from Morse theory and the analysis of asymptotic 
fluctuations of critical points of random eigenfunctions 7j. The expressions (1.4)-(1.7) are supported by 
extensive numerics [Sj. 

Remark 2. Building upon previous results [l9j . we are now able to present a full characterisation for the 
asymptotic behaviour for the variance of the three Lipschitz-Killing curvatures £j, * = 0,1, 2, for the excursion 
sets of random spherical eigenfunctions on S 2 . In this setting these three LKC’s correspond, respectively, to 
the EPC (£o)> half the length of level curves (£i), and the excursion area (£ 2 )- Indeed, it was shown [13] 
that the variance of the excursion area for spherical Gaussian eigenfunctions satisfies 


lim £ Var[£ 2 (A u (/ £ ; S 2 ))] = u 2 cf 2 {u) = [H^u) 
£—>•00 


H'iutfHuf. 


Note that the asymptotic covariance in (1.41 can be positive, negative or null depending on the choice of the 
intervals I\ and 1 2 (see Figure [l]). From (1.41 and (1.5) it follows also that, for every intervals Ii, J 2 C R 
such that the corresponding variances do not vanish, as I goes to infinity, %(A / 1 (fe ; <S 2 )) and y(A / 2 (/£; S 2 )) 
are asymptotically perfectly (positively or negatively) correlated, i.e. 

Corollary 2 . For all intervals /1, J2 such that 

Var[ X (A /l (/,;5 2 ))],Var[ X (A / 2 (/,;5 2 ))] ^ 0, 

as £ —> 00, 

|Corr[x(A/ 1 (ff, S 2 )) , X (A /2 (/,;5 2 ))]| = 1 + O (r 1 / 2 ). 

A similar form of degeneracy was earlier observed for level curves in m- From Theorem |T] we also have 
the following corollary for half-intervals Ji = [ui,oo) and I 2 = [112,00) (see also Figure [ 2 ] and Figure [3]): 

Corollary 3 . As £ —> 00, for U\,U2 £ R, 

( 1 . 6 ) 

In particular, if u\ = U2 = u, we can present an analytic expression for the variance: 

(1.7) Var {x(A u (ff,S 2 ))} = H 3 (u ) + 2 H 1 (u)} 2 e ~ u2 + 0(£ 5 / 2 ), 

where H q {-) are the Hermite polynomial of order q. 


Cov[x(A Ul (fr,S 2 )), x(A„ 2 (/^;5 2 ))] = ^uiu 2 (u? - l){u\ - l)e ?e * +0(£ 5/2 ). 
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(a) positive range 


(b) negative range 


u 2 u 2 

Figure 2. ^uiu 2 (ul — 1 )(u| — l)e - ^e - ^ 



Figure 3. ^[H 3 (u) + 2H 1 (u)] 2 e~ v 


On the other hand !3Sj formula (18) (see also [IS]) asserts (in a slightly different form) that, for the variance 
of the boundary length of excursion sets, the following result holds 

(1.9) lim f -1 Var [Ci(A u (ff 1 S 2 ))] = const x i t 4 </> 2 (u) = const x [H 2 (u) + H[(u)] 2 ft 2 (u). 

£.—^-oo 

Likewise the asymptotic variance of the Euler-Poincare characteristic, derived in Corollary [3j may be written 
as 


lim i 3 Var[£ 0 (A u (/^;5 2 ))] = hw 3 - w] 2 c/> 2 (u) = \[H 3 {u) + H' 2 {u)] 2 (j) 2 (u). 
£—>oo 4 4 


( 1 . 10 ) 

We may unite the asymptotic expressions for the variance of the first three Lipschitz-Killing curvatures in 


(1.81, (1.9) and (1.10) into a single formula: 


( 1 . 11 ) 


lim e 2k 3 x Var [C k (A u (ff, S 2 )\ = const x [H 3 _ 2k {u) + H 2 _ 2k (u)} 2 (t> 2 (u), k = 0,1,2. 




A comparison of (1.11) with expressions 0, ( |1.2[ ) and ( |1.3[ ) seems to suggest the existence of an (as¬ 
ymptotic) second order Gaussian Kinematic Formula for spherical Gaussian eigenfunctions. We leave the 
investigation of the general validity of such en expression for higher dimensional spheres to future research. 
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Remark 3. For all zz £ R, we have 

x(A u (fc,S 2 )) / M 

E[x(A u (fr,sm 

with the usual convention X n = O p (a n ) meaning that the sequence \X n \/a n is bounded in probability; i.e., 
in the high frequency limit £ —> oo, the ratio of the realised and expected value for the EPC of the excursion 
will converge to unity in probability for all u £ R. 


2. Background on Morse theorem and (approximate) Kac-Rice formula 

2.1. Morse theorem. We start by recalling a general expression for the EPC by means of so-called Morse 
Theorem (see P Section 9.3). Assuming that Ad is a C 2 manifold without boundary in and that 
h £ C 2 ( Ad) is a Morse function on Ad (i.e. its Hessian is non degenerate at the critical points), we have 


( 2 . 1 ) 


dim(A-f) 

x(M)= (-^ViCAd,/!), 

3= 0 


where /j,j( Ad, h ) is the number of critical points of h with Morse index j, i.e., the Hessian of h has j negative 
eigenvalues. In order to develop our results we will need to exploit (2.1) in the case of excursion sets of 
spherical eigenfunctions; to this end, let us first recall some basic differential geometry on S 2 . The metric 
tensor on the tangent plane T(S 2 ) is given by 


g{9,v) = 


0 

sin 2 


For x = {6,<p) £ S 2 \ {N, 5} (N,S are the north and south poles i.e. 9 = 0 and 9 = tt respectively), the 
vectors 

_ <9 _ 1 <9 

68 39 ’ Cip sin 9 dtp ’ 

constitute an orthonormal basis for T X (S 2 ); in these coordinates the gradient is given by V = (J^, 

The Hessian of a function / £ C 2 (S 2 ) is the bilinear symmetric map from C 1 (T(S 2 )) x C 1 (T(5 2 )) to C 0 ^ 2 ) 
defined by 

V 2 /( X, Y ) = XYf - V x Yf, X, Y £ T(S 2 ), 

where Vx denotes Levi-Civita connection (see e.g. [2| Chapter 7 for more discussion and details). For our 
computations to follow we shall need the matrix-valued process X 2 E ft{x) with elements given by 

E f ^(**0} a, b=G = {(^ /^(*^))(^ai ) }a.&— 6,(p i 

where E = {eg,e v }. In coordinates as above, this matrix can be expressed as 

jY _ r 0 _d_ _pV 3 _d_ 1 r d 2 

dO 2 06 dO 06 dip sin 6 L < 

r d 2 


V^(*) = 


_ _ _p^ d _ pO d 1 

sin 0^-dOdip ipO dip 1 OipdOi 


d 2 
dO 2 

1 r d 2 _ cos 0 d l 

sin O'- dOdip sin# dtp J 


_ pV _d _ pO d ] 

sin 6 dOdip <pO dip Oip dO -I 

1 r d 2 pep d pO d ] 

sin 2 6 L dip 2 ^ W dip 1 W dO J 


d z 


cos 6 d 
sin 6 dtp J 


sin# L dOdip 

s!iAe[^ + sin0cos0 MJ 


Here p b are the usual Christoffel symbols, see e.g. [5] Section 1.1, which allow to compute the Levi-Civita 
connection: 


^£ a e b — r ab e S + rfey, 

More explicitly, Christoffel symbols for S 2 are given by 


a,b = 6, ip. 


-pC7 _ ptf 

1 Oip — 1 06 


— p<P — pV — n 

— 1 ipip — L eo ~ 


Y 6 = — sin 0 cos 0 , 


T* g = cott 


We now state the Morse representation for the Euler characteristic of the excursion set: let Ad and h in 
( |TT| ) be A/(/ f ;5 2 ) and f e \ Al (/f; 5 2 ) respectively, we have 


( 2 . 2 ) 


X(MfcS 2 )) = £(-!)>,., 


3=0 


where 


lij = #{x £ S 2 : f e (x) £ I,Xf e (x) = 0, Ind(-V|/ £ (x)) = j}, 
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Ind(.M) denoting the number of negative eigenvalues of a square matrix M. More specifically, g o is the num¬ 
ber of maxima, g i the number of saddles, and g 2 the number of minima in the excursion region Ai(fg,S 2 ). 

2.2. Kac-Rice formula. The Kac-Rice formula is a standard tool (or meta-theorem) for expressing the 
(factorial) moments of the zero crossings number of a Gaussian process in terms of certain explicit integrals. 
In our case, we are interested in counting the critical points of fi, i.e. the zeros of the map x —>■ V fe(x). Let 
£ C R n be a nice Euclidean domain, and g : £ —> R n a centred Gaussian random field, a.s. smooth. Define 
the 2 -point correlation function of critical points Iv 2 = /\ 2; g ■ £ 2 —> R 

K 2 (x,y) = <P(v g (x),v g ( y )){ 0,0) • E[|detV 2 g(ir)| • \det\/ 2 g(y)\\\/g(x ) = V g{y) = 0], 

where <^(Vg(x),Vg(y)) the Gaussian probability density of (X7g(x),X7g(y)) £ R 2n . Let J\f c (g) = J\f c (g,£) = 
#{x £ £ : S7g(x) = 0}; by virtue of [S’ Theorem 6.3, we have 


E[M c (g,£)-(U c ( 9l £ )-!)] = 


£x£ 


K 2 {x,y)dxdy, 


provided that the Gaussian distribution of (V g(x),Vg(y)) £ R 2 " is non-degenerate for all (x,y) £ £ 2 , on 
the validity condition of Kac-Rice formula in the Gaussian case, see [5) Theorem 6.3 and Proposition 1.2, 
and [T6] Section 1.4. Moreover for V i,7 ? 2 C £ two nice disjoint domains, under the same non-degeneracy 
assumptions for every (x,y) £T>i x P 2 , we have 


E[A^ c ( ff ,2?i) -Ar c (g,V 2 )] = [[ K 2 (x,y)dxdy. 

J J X^i x Do. 


Jv ixX >2 

It is easy to adapt the definition of the 2-point correlation function in order to investigate, for example, the 
maxima with values lying in an interval ICR: we re-define /\ 2 as 

-^2,0,0 (x,y,I,I) = </ ) (Vg(x),Vg(y))(0'0) 

x E[|detv 2 5 (x)| • |detV 2 0 (?/)| ■ 1 7 ( 3 ( 0 ;)) • 1 i(g{y)) ■ l{ind(-v 2 g(x))=o} ■ l{ind(-v=g(g))=o}IVg(x) = Vg{y) = 0], 

where 1 / is the characteristic function of I on R. 

For the Kac-Rice formula on manifolds we refer to j2] Theorem 12.1.1, in particular, let 

G S 2 : fi{x) £ /, V//(a:) = 0} 

be the total number of critical points in / of {fe(x),x £ <S 2 }; we have 

■ (A/?(/*) - 1)] = JJ K 24 {x, y, I, I)dxdy : 


(2.3) 

where 


S 2 xS 2 


(2.4) 


K 2 ,e{x,y,I,I) = <j>(Vfe(x),Vft(y))( °> 0 ) 

x E[| detV 2 E fe(x)\ • | det V|/^(y)| ■ 1 /(/?(*)) • h(fe(y))\V fa(x) = V f e {y) = 0], 

One technical difficulty in working with the spherical Gaussian eigenfunctions fi in ( |2.3| ) is related to the fact 
that the Gaussian distribution of (fe{x),Vfe(x),V%fg(x)) is always degenerate. However, this issue can be 
handled by writing fe as a linear combination of second order derivatives, and thus reducing the dimension 
of the Gaussian vector involved in the evaluation of K 2 , see [7]. 

A much trickier issue arises when we need to validate a sufficient non-degeneracy assumptions due to the 
technical difficulties of dealing with 10 x 10 matrices depending on both x and y (and £). Following [7] and 
m, we do not claim the (precise) Kac-Rice formula ( |2.3| ) but rather an approximate version, see [7] formula 
(3.5), equivalent to (2.3) up to an admissible error. 

First note that, by isotropy, K 2 {x,y) = AT 2 (d(x, y)) depends only on the (spherical) distance d(x,y) = 
a,rccos((x,y)) between x and y. In view of this, we note that it is convenient to perform our computations 
along a specific geodesic; in particular, we constrain ourselves to the equatorial line 9 X = 9 y = 7 r/ 2 ; it is 
immediate to see that here the gradient and the Hessian are 

■ — t d ® \ 

\e=^/ 2 ~ (gff’fy)’ 



r a 2 

ae 2 

a 2 1 

dOdp 

1 6»=tt/2 

a 2 

a 2 

_ dQdtp 

dip 2 


The basic idea is to split the range of integration in (2.3) into two parts: the “short range” regime 
d(x,y) < C/£ and the “long range” regime d(x,y) > C/i, C denoting a sufficiently big positive constant. In 
the short range regime Kac-Rice formula holds only approximately, but, by a partitioning argument inspired 
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from m (see also iTSj), it is possible to prove that its contribution is 0{£ 2 ). In the long range regime 
d(x,y ) > C/£ the Kac-Rice formula is precise. The above yields 

(2.5) — 1)] = f ff K 2/ (x,y;t 1 ,t 2 )dt 1 dt 2 dxdy + 0(£ 2 ), 

Jd(x,y)>CIi J Jlxl 

where 

(2.6) K 2 ,e(x,y;ti,t 2 ) 

= Tx,yAh,t 2 , 0 , 0 )IE[| detV|/ f (a;)| • \ det V 2 E f £ (y)\\Vfi(x) = V/ f (y) = 0 ,f t (x) = t ly f e (y) = t 2 ], 

and ip x ,y,i is the density of the 6 -dimensional vector (fe(x), fi(y), V fi(x), V/^(?/)). For further details on the 
proof of (2.5), see [7| Section 3.4.1 and Section 3.4.2. 

As it will become clear from the proof of Proposition [2] below, we obtain a considerable simplification 
in our calculations since, during the application of the (approximate) Kac-Rice formula for studying the 
variance of the EPC, we can get rid of the absolute values in ( |2.6[ ); in fact, for g as before a smooth, centred 
Gaussian random field, we observe that (0 Lemma 4.2.2) 

(-1) J | det V| ; 5 , (a:)|]l {Ind( _ V | 9 ( a . )) _ J -} = (-1) J sgn(det V 2 E g{x)){det V 2 E g(x))l {lnd{ _ V 2 Eg{x))=j} 

det( V E g[x) ) l{Ind( — V^.g(x))=j} > 

since sgn(detV|g(a;))l {Ind( _ V | g(a . ))=j} = (-1 ) J+1 and -det\7 2 E g(x) = det(-V| g(x)). Hence 

2 

(2.7) l det V|c/(a:)|l{ind(-v| ;S (x))=n = det(-V|ff(a:)). 


7=0 


3. Proof of Theorem [T] 


Let /» CM, i = 1,2 be two interval in the real line; in the argument to follow we shall adopt the following 
notation: 


e 52 : M x ) G = 0 ,lnd(—V| fe(x)) = j }, 


3 = 0,1,2, i = 1,2. 


Theorem [Tj is a straightforward application of Proposition [T] and Proposition [2j The first building block is 
the approximate Kac-Rice formula for covariance computation: 


ld(x,y)>C/l J J IiXl 2 


J 2 ,t(x, y; ti,t 2 )dtidt 2 dxdy + 0(£ 2 ) 


Proposition 1. There exists a constant C > 0 sufficiently big, such that 

( 3 . 1 ) £ (-1 y +k nn r Af e )y k ; 2 (h)} = f 

j,k —0 Jd ^ 

where 

( 3 . 2 ) 

J 2 ,i(x,y;t ly t 2 ) 

= Tx,y,e{ti,t 2 ,0,0) E[det(-V|/e(x)) • det(-V|/e(y))|V/e(a:) = Vfa(y) = 0 , f e (x) = ti,fe(y) = t 2 \. 

Proof. We start by observing that 

2 2 

( 3 . 3 ) £ (-iy +k E\^{hffi k , 2 {fi)\ = Y.(~ l ) j+k ^rAh)»M)] 

j,k =0 j^k 7=0 


For the non diagonal terms in (3.3) with j ^ k, j, k = 0,1, 2, we directly obtain (see [7] Section 3.4) that for 
any sufficiently big constant C > 0, we have 

E [/H;i(/e)/bfc;2(/e)] = / // K 2 ,t l j,k{x,y;t 1 ,t2)dt 1 dt 2 dxdy + 0 (£ 2 ), 

Jd(x,y)>C/£ J JI 1 XI 2 

where 

K2,ej.k{x,y;ti,t 2 ) = <p x ,y,e(ti,t 2 ,0,0) 

E[| det V 2 E f((x) \ ■ |detV|/f(y)| • l {Ind( - V |/,(x )) =j } ' l{ind(-v|/R 3 /))=fc} | Vfe(x) = V f e (y) = 0 , f e (x) = t u f t {y) 
To work out the diagonal terms j = 0,1, 2, in ( |3.3| ), we introduce the following notation: 

Vj-Ai'ife) = #{x G S 2 : fe(x) £/,\ h',Vf e {x) = 0,Ind(-V| f t {x)) =j}, 
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9 j-,ini'(ft) = #{x e S 2 : f{(x) e I,n Ii',Vft(x) = 0,Ind(-V|/*(a;)) = j}, 
with i, %' = 1, 2, so that 

f J ’j]l(ft)fJ'j-,2(ft) = (fJ>j-l\2(ft) + Mj;ln2(/e))(Mj;2\l(/e) +Mj';ln2(/^)) 

and then 

= E[/i J - ; i\ 2 (/6)M i j;2\i(/«)] + ®'[Mj;i\2 (^ig)Mj;in2 (/e)] + E[/ij;2\i(/«)Mj;in2(/e)] + ^■[V’j-,in2(fe)V'j-,2ni(fe)} 

= ^[f2j-l\2(ft)^j-2\l(ft)] + ^[Pj-,\\2(ft)l2j',ln2(ft)] + ;2\l( ft) ;ln2( ft)\ 

(3.4) + ^[i2j ; \ri2(ft)(l2j;in2(fe) — 1)] + ^[Pj;in2(ft)]- 

For the last term in (3.4) we note that the expected value of the EPC of the excursion set is 0(£ 2 ), while 
for the other terms we can apply again the approximate Kac-Rice formula. For example we have: 

K 2 ,e,j,j(x, y; t 1 ,t 2 )dt 1 dt 2 dxdy + 0(£ 2 ), 


E[% ; in2(/e)(M.j;in2(/e) — 1)] = 


1 d(x,y)>C/t J Jh\I 2 Xl 2 \h 


'd(x,y)>C/£ J J 7in/2X/in/2 


K2,tj,j(x,y,t 1 ) t 2 )dt 1 dt 2 dxdy + 0 (£ 2 ), 


and then 


E [N-Aft)VrMt)\ = 


Id(x,y)>C/t J J It Xl 2 


^2,t,j,j(x, y, ti 1 t 2 )dtidt 2 dxdy + 0 (£ 2 ). 


d(x,y)>C / 1 J JI\Xl 2 


/ / J2,t(x, y; ti,t 2 )dtidt 2 dxdy + 0(£ 2 
J J It X In 


We can apply now the identity (2.7) to get: 

2 

E (-iy +k n^i(h)^ 2 (h)} = 

j,k =0 

where 

J 2 ,t(x,y;t i,f 2 ) 

= Vx,y,t(h,h,0,0) E[det(-V| f e (x)) • det(-V % fe(y))\Vfe(x) = Vfi(y) = 0 ,fc(x) = h,f e (y) = t 2 }. 

Our second tool yields an analytic expression for the alternating sum in the variance computation. 

Proposition 2. 

2 

E (-l) J+fc E[^ ;1 (/,) Mfc;2 (/,)]-E[ X (A /l (/,;5 2 ))]E[ X (2l /2 (/,;5 2 ))] 

j,k=0 

£ 3 r /' f 

Pi(ti)dti / Pi(t 2 )dt 2 - 


□ 


4 L 


h 


>I 2 


92 (ti,t 2 )dt 1 dt 2 + 16 / g3(h)dti g 3 (0)dt 2 +0(£ 5/2 ), 


hxh 


' h 


Hi 


where 


(3.5) pi(f) = 


* / (xitVS- x\ - a:^) exp 1 -\t 2 1 exp{-ha; 2 + x\ - \Z8tx 1 )}dx 1 dx 2 

(2tt) 3 /^ J r 2 l 2 J 2 


(3.6) 9 2{t 1 ,t2) = ~ ) , 


L / 'Ll ( ziV ~ 8ti zl ) exp D exp {L + Z2 

x (wiV 8 t 2 - w\ - w'^j exp j-^t 2 j ex P + w\ - %/ 8 t 2 wi)| 

-6 + (3*i - V^zi) 2 + (3f 2 - V2w\) 2 


dzidz 2 dw\dw 2 , 


and 

(3.7) 


93(t) = \ j^ 3/2 J r2 (ziVst - zl - z£) ex P |-y 2 |exp|-i(z? + z% - v^tzi) j 3 - (3i - v^zi) 5 


dzidz2- 
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Proof. In view of Proposition [l] and by isotropy we have to study the asymptotic behaviour of 

fir/2 


(3.8) 


107T 


>c/e 


sincf dcf -E[x(A Il (f e ;S 2 ))]E[x(A l2 (f l ;S 2 ))} + 0(r). 


hxi 2 


Again we stress that J 2 ^ in (3.2) is analogous to A ' 2 1 in (2.4) except for the fact that the absolute value of 
the Hessian determinant has been dropped (by means of Morse theorem). 

The proof of this proposition follows along the same lines as in the argument given in 7] Section 4.1.2 
where we study the asymptotic behaviour of 

16n 2 [ [[ K 2 /{(t>-,t ll t 2 )dt 1 dt 2 sin(j) (E [A / - /(/*)]) 2 

Jc/e JJixi 


to obtain the variance of the number of critical points. Therefore here we just sketch the main steps and we 
refer to [7] Section 4.1.2 for a complete proof. 

The asymptotic analysis is based on the properties of multivariate conditional Gaussian variables, and 
on an asymptotic study of the tail decay of Legendre polynomials and their derivatives that appear in the 
conditional covariance matrix of the Gaussian vector. In fact, for d(x,y) > C/£, C large enough, Kac-Rice 
formula holds exactly and we one can exploit the fact that a Gaussian expectation is an analytic function 
with respect to the parameters of the corresponding covariance matrix outside its singularities. It is then 
possible to compute the Taylor expansion of these expected values around the origin with respect to the 
vanishing entries 


a = a t {(j)) = ( 01 ^( 0 ), a 2 , t {<t>), a 3 ,e(4>), ov(0)> a 6 , *(</>), a 7< e (</>), a s ,e{^)) 


of the conditional covariance matrix Af(</>) = A(a) (see [7] Appendix B) of the centred Gaussian random 
vector 

^(V 2 //(*), V 2 /,(y)|V/ f (x) = Vf e (y) = 0). 

Three terms in the Taylor expansion (depending on the intervals I\ and I 2 ) give an asymptotically significant 
contribution, whereas the rest is negligible: 


I6n 2 [ [[ J 2/ {(t)]ti,t 2 )dt 1 dt 2 sin(j)d(j)~E[x{A Il {ff 1 S 2 ))\E[x(A l2 (ff 1 S 2 ))] 

Jc/i JJhxIo. 

14 


| / Pi{ti)dt! I pi(t 2 )dt 2 — 16 


(3.9) 


32 


h JI 2 

A, (a; (l , t2 ) 


hxl 2 


hxl 2 

dtidi 2 } + 0(f 5 / 2 ), 


d , 

— q{a;tut 2 ) 

lda 3 


dt\dt 2 


here we set 


g(a;ti,t 2 ) 


= TWW f f (ziV8t 1 -z\- zfj ■ (w 3 v / 8t 2 -w\- w%) 

\ 2tt) j jjk 2 xR 2 V /V / 

x q( a; t\, t 2 ; z 1; z 2 , w 1 ,w 2 )dz i dz 2 dw 1 dw 2 , 


with 


q(a;t 1 ,t 2 -,z 1 ,z 2 ,u>i,w 2 ) = 


exp \ —\v tx ,t 2 z 2 , wi,w 2 )A.(a) 1 v tlM {zi,z 2l w 1 ,w 2 ) i 


v/det(A(a)) 


and pi defined in (3.5). Note that the zeroth order term in the Taylor expansion cancels out with 

E [ X (A Jl (/,;5 2 ))]E[ X (A 72 (/,; < S 2 ))] 

that is of order 0(£ 4 ). The expressions for g 2 and g 3 in ( |3.6[ ) and ( |3.7| ) follow from the evaluation of the 

□ 


partial derivatives in formula ( |3.9| ); once more we refer to [7j Section 4.1.2 for details. 
We can now prove Theorem [T] 


Proof of Theorem^ 7} We first write: 

Cov[x(A / 1 (f e ;S 2 )), X (A l2 (f t ; S 2 ))] = E[ X (A h (/,; S 2 )) X (A h (/*; S 2 ))} - E[ X (A 7l (/,; <S 2 ))]E[ x (A /2 (/,; S 2 ))] 
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where, in view of (2.2), 

2 2 

EM4(//;5 2 )W4(/i;5 2 ))]=E[ E (-i) j+ W(/*)a*;2(//)] = E 

j,k =0 j,k —0 

Now by Proposition [2] the covariance is asymptotic to 

Cov[x(A 7l (fe ; 5 2 ))x(^4/ 2 (fi ; S 2 ))] 

f3 r f /• 


4 L J 

Now define 


p 1 {ti)dt 1 / pi(t 2 )dt 2 - 


g 2 (t 1 ,t 2 )dt 1 dt 2 + 16 / P 3 (ii)rfii / g 3 (t 2 )dt 2 + 0(£ 5/2 ). 


11 X/o 


p 2(0 = ^ 2 ^ 3/2 J (3t — v / 2x 1 ) 2 (x 1 t\/8 — x\ -£ 2 )exp j-^t 2 jexp j-*(:r 2 + a;2 - j dx 3 dx 2 , 

it is easy to see that the functions g 2 and g 3 in (3.6) and (3.7) can be rewritten as 

g2(tl,t 2 ) = —3pi (tl)pi(t 2 ) + ^P2(tl)pi(t 2 ) + ^Pl{tl)p 2 (t 2 ), g 3 (t) = jjpi(t) - ^P2{t). 

Moreover p 3 and p 2 can be explicitly computed and we have: 


. \ a/2 , _tl 

Pi(t) = —j= (t 2 - l)e 2 , 
\h r 


P2(t) = ^L(t 4 +t 2 - 4)e ' 2 . 

V7T 


It follows that we can rewrite the coefficient of the leading term in the following form: 

[ p 1 (t 1 )dt 1 [ pi(t 2 )dt 2 - [[ g 2 (t 1 ,t 2 )dt 1 dt 2 + 16 f g 3 {t{)dt 3 f g 3 (t 2 )dt 2 

Jh Jl 2 JJhXl2 dll ■' r2 

= I li il 2 ,i - [[ [—3pi(ti)pi(t 2 ) + \p2{t\)p\{t 2 ) + ^.Pi{ti)p 2 (t 2 )\dt 1 dt 2 

Jjiixi 2 1 l 

+ 16^ [§Pi(*i) - \p2{ti)]dt 3 J 1 l^pi(t 2 ) - i p 2 (t 2 )}dt 2 
= x 1 , 1 x 2 , 1 - [-3 X m x 2i1 + ^X lj2 x 2>1 + ^X M X 2 , 2 ] + le^x^ - |x 1;2 ][|x 2>1 - ^x 2;2 ] 


- 25 x X 5 x X 5 r X +E X 

— _ ^ _ -^l,l - ^2,l — ~Xi ? 2-^2,1 — ~Xi j iX2 J 2 i ~X]_ ? 2-^2,2 > 


(3.10) 

where Xjj, for *, j = 1,2, are given by 
Xy J - I Pj(t)dt 


. , a/ 2, 9 t 2 

Pi(t)=-^(t 2 -l)e 2 , 

\/7T 


p 2 {t) = ^(t 4 + t 2 - 4)e” 

\/7T 


Formula (3.10) can be further simplified as follows: 

1 


— / +Ati-\)e-*dh / (-^ + 4t 2 -l)e-^dt 2 . 


27T 


'/l 


□ 


Proof of Corollary [TJ The asymptotic expression for the variance in (1.5) follows immediately by setting 
Ii = I 2 = 7, and we have: 


V2n ■ 


(-1 4 + 4t 2 - l)e“^di|“ + 0(f 5/2 ). 


Var[x(A/(/^;5 2 ))] = j 

That is the statement of Theorem E □ 

Proof of Corollary^ In the particular case where I\ = [iti,oo) and I 2 = [it 2 ,oo), we have the following 
explicit form for the leading term of the covariance: 

u 2 u 2 

X\I 2 = u\u 2 {u\ — l)(it 2 — l)e - ^ . 

Also, for Ii = I 2 = [u, 00 ), our expression reduces to 

Var [ x (A u (ff,S 2 ))\ = ^(u - u 3 ) 2 e~ u2 + O^ 2 ) = £-(H 3 (u) + 2H 1 (u)) 2 e~ u2 + 0(£ 5/2 ), 


87 T 


87T 


as claimed. 


□ 
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